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We study the electron motion near magnetic field steps at 
which the strength and/or sign of the magnetic field changes. 
The energy spectrum for such systems is found and the elec- 
tron states (bound and scattered) are compared with their 
corresponding classical paths. Several classical properties as 
the velocity parallel to the edge, the oscillation frequency per- 
pendicular to the edge and the extent of the states are com- 
pared with their quantum mechanical counterpart. A class of 
magnetic edge states is found which do not have a classical 
counterpart. 

73.40-c; 73.50-k; 73.23 



I. INTRODUCTION 

The transport properties of a two-dimensional electron 
gas (2DEG) subjected to a nonhomogeneous perpendic- 
ular magnetic field (periodically modulated or not) have 
been the focus of a great deal of research in recent years. El 
Current fabrication technologies permit to create nonho- 
mogeneous magnetic fields on a nanometer scale by delib- 
erately shaping or eiujving the 2DEGJ3 or by integration 
of superconductingem or ferromagnetic materialsoB on 
top of the 2DEG. This will add a new functional dimen- 
sion to the present semiconductor technology and will 
open avenues for new physics and possible applications.Q 

Theoretically, the effect of nonhomogeneous magnetic 
fields on a 2DEG have been studied both in the ballistic 
and the diffusive regime. The resulting perpendicular 
magnetic field can act.-as a scattering centre Jju'u but can 
also bind electrons ,0113 and so influence the transport 
properties of the 2DEG. In transport calculations one 
needs the electron states, which are obtained by solving 
the Schrodinger equation. 

MulleitJ studied theoretically the single particle elec- 
tron states of a 2DEG in a wide quantum waveguide un- 
der the application of a nonuniform magnetic field and 
showed that in the case of a magnetic field modulation in 
one direction, transport properties also become one di- 
mensional and electron states propagate perpendicularly 
to the field gradient. 

Making use of this decoupling, the electron states for 
different nonhomogeneous magnetic field profiles along 
one dimension were investigated, i.e. for a periodically 
modulated magnetic fieldlinlj, for magnetic xuiaiitum 
steps, barriers and wells in an infinite 2DEGLtE3 and 
in a narrow waveguide. E3 



In this paper we consider an infinite 2DEG subjected to 
a step-like magnetic field, i.e abruptly changing in mag- 
nitude or polarity at x = 0, in one dimension (taken to 
be the ^-direction). Preliminary results were presented 
in Ref. ^0|. First the situation for two opposite homo- 
geneous magnetic fields with the same strength will be 
considered. The classical trajectories correspond to snake 
orbits and were already used in the seventiesEU to describe 
electron propagation parallel to the boundary between 
two magnetic domains. Back then, one was interested 
in understanding the electron transport through multi- 
domain ferromagnets and it turned out to be more con- 
venient to work with the classical trajectrories than with 
the corresponding electron states, which allows one to use 
a semi-classical theory which reduces the complexity of 
the theory considerably. We are interested in transport 
through a 2DEG situated in a semiconductor in which 
the Fermi energy is orders of magnitude smaller than in 
the metallic systems of Ref. 5jl 

We will study thoroughly the quantum mechanics of 
such electron states in a 2DEG subjected to this step 
magnetic field profile, and we will compare them with 
their classical counterpart. We will discuss the energy 
spectrum and the corresponding electron states, and de- 
rive several properties. We will show the existence of 
states wich have a velocity opposite to the expected clas- 
sical orbits. Additionally, we will show that adding a 
background magnetic field modifies the spectrum and the 
states considerably. 

The paper is organized as follows. In Sec. II we present 
our theoretical approach. In Sec. Ill we calculate the en- 
ergy spectrum, the wavefunctions and their correspond- 
ing group velocity, and compare this with their quantum 
mechanical counterpart. In Sec. IV we study the influence 
of a background magnetic field on the quantum mechan- 
ical and classical behaviour. In Sec. V we focus on the 
negative velocity state, and finally, in Sec. VI, we con- 
struct time dependent states, and interpret them classi- 
cally for several magnetic field profiles. 



II. THEORETICAL APPROACH 

We consider a system of noninteracting electrons mov- 
ing in the xy-planc in the absence of any electric 
potentials. The electrons are subjected to a magnetic field 
profile B = (0,0, B z (x)). First, we will study the elec- 
tronic states near the edge of two magnetic fields with 
opposite strength 
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B z (x) = B [26(x) - 1] 



(1) 



which is independent of the y-coordinate. Next, we will 
consider the influence of a background magnetic field B 
on these states, which results in the magnetic field profile 



B z (x) = B [26(x) 



l]+B. 



(2) 



In the following we will use B l = B z (x < 0) and B r = 
B z (x > 0), to denote respectively the magnetic field on 
the left and the right hand side of the magnetic edge. 

The one-particle states in such a 2DEG are described 
by the Hamiltonian 
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Taking the vector potential in the Landau gauge, 

A = (0,xB z (x),0), 
we arrive at the following 2D Schrodinger equation 



(3) 



(4) 



dx 2 
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d_ 

dy 



+ ixB z (x) 



where the magnetic field is expressed in B , all lengths 
are measured in the magnetic length Is — \/Tic/eBo, 
energy is measured in units of Tiu Cl with uj c = eBo/m e c 
the cyclotron frequency and the velocity is expressed in 
units of Ib^c- H and p y commute due to the special form 
of the gauge, and consequently the wavefuntion becomes 



ip(x,y) 



- ikv <f>n,k{x), 



(5) 



which reduces the problem to the solution of the ID 
Schrodinger equation 



2dx 2 



V k (x) 



4>n,k(x) = E n<k cj) nik (x), 



where it is the fc-dependent effective potential 
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Vfc (x) = - [xB z (x) - ky 



(6) 



(7) 



which contains the two dimensionality of the problem. We 
will solve Eq. (Q) numerically by use of a discretization 
procedure. In some limiting cases analytical results can 
be obtained. 



III. IN THE ABSENCE OF A BACKGROUND 
MAGNETIC FIELD 

Let us first consider the case when no background mag- 
netic field is present. The situation is then symmetric, 
and more easily to solve. The effective potential for this 
case is shown in Fig. ||(a) for k = — 2 (dotted curve) and 



k = 2 (solid curve) . We notice from Eq. (Q) that this po- 
tential is built from two parabolas, with minima situated 
at x l = —k, and x r = k, thus respectively on the left and 
right hand side of the magnetic edge. The total potential 
has for k > two local minima respectively at x = — k 
and x — +fc, while for k < it has only one minimum 
at x = 0. Before we describe the energy spectrum of the 
snake orbits and their corresponding properties, we first 
discuss the limiting behaviour. 



A. Limiting behaviour for k — > ±oo 

For k — > oo, the minima of the parabolas are sit- 
uated far from each other. The electrons are in the 
Landau states of two opposite magnetic fields, one on 
the left, the other on the right, and they are not in- 
teracting with each other. The electron wavefunctions 
are given by (L\x) — C rn H rn (x + k)e~( x+kS> I 2 , and 
(R\x) = C m H m (x — fc)e~( x ~ fc ) / 2 , respectively, where 
H m 0*0 is the Hermite polynomial. For decreasing k the 
parabolas shift towards each other, and the electrons 
will start to "feel" each other. In terms of wavefunc- 
tions, this results in a parabolic cylinder function 4>(x) = 
D E _i [v2 (x — fc)] matched at x — 0, with the condition 



that 4-Dz 



(«)U- 



V2k 



or D E _ 1 {-V2k) 



0, 



for the symmetric and the antisymmetric wavefunction, 
respectively. This leads to a change in energy of the 
electron states, which can be understood as a lifting 
of the degeneracy of the two original electron wave- 
functions. The energy can then be written as E±(k) — 
(L\ H \L) (k) ± (L\ H \R) (jfe) = E (jfe) ± AE (k) with the 
corresponding wavefunctions \<f>) = \R) =p \L). One can 
see that the presence of the second parabola results in 
two effects: (1) a decrease of E (k) due to the finite pres- 
ence of the wave function in the other parabola, and (2) 
a splitting of the energy level due to the overlap, i.e. one 
level (E + ) shifts upwards, while the other {EJ) shifts 
down. For k — > oo, this results in the following first-order 
approximation to E and AE: 



E m (k) = i + m - 



'2m-l-k 2 



AE m (k) 



>2m+l 



=e-« k 



(8a) 
(8b) 



In the other limit k — > — oo, the effective poten- 
tial can be approximated by a triangular well V {x) — 
k 2 /2 — kB x. Solutions for this potential consist of Airy 
functions, again matched at x — with the condition 
that (/>' (0) = or (0) = which results, respec- 
tively into the anti-symmetric wavefunction 4>2 m (%) — 



C 2m (k) (\x\ / x) Ai ZAi', m +i + (2k) 
metric 

one (j) 2m +i (x) = C 2m +\ (k) Ai 
respectively with energy 
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E 2m (k -> -oo) = - k 2 - z A i',m+i (2 \k\) 
E 2m +i (k -> -oo) = - A; 2 - ZAi, m +\ (2 |fc|) 



2/3 
2/3 



(9b) 



where z A%n (= -2.338, -4.088, 

-5.521, - [3tt (4n - 1) /8] 2/3 ) and z^i'.n (= -1.019, 

-3.248, -4, 820, - [3tt (An - 3) /8] 2/3 ), denote respec- 
tively the n th (n = 1, 2, 3, oo) zero of the Airy function 
and of its derivative. One can see that for increasing neg- 
ative k, the difference between the two energy branches 
increases, which is to first order linear in \k\. Namely the 
more negative k, the narrower the well, thus the more 
the energy levels are shifted from each other. 



B. Spectrum and velocity 

Solving Eq. (^) numerically gives rise to the energy 
spectrum shown (solid curves) in Fig. ||(a). For k = oo, 
we obtain the earlier mentioned Landau levels, which 
are labelled with the quantum number m. Each level is 
twofold degenerate. For decreasing k, the degeneracy is 
lifted and they separate into two different branches with 
eigenstates |2m) and |2m + 1) and eigenvalues Ei m and 
£2m+i, and corresponding quantum numbers n = 2m 
and n = 2m + 1. This quantum number n does not only 
result from arranging the levels according to their lowest 
energy, starting with n — 0, but it also reflects the num- 
ber of nodes of the corresponding wavefunction. Notice 
that the levels have now a non-zero derivative, i.e. elec- 
trons propage in the y-direction, and their group veloc- 
ity is given by v (k) = —dE (k) jdk. (The minus sign 
appears here because in Eq. (||) we took k y = —k). 
This group velocity is plotted (solid curves) in Fig. ||(a) 
for the 6 lowest levels. For k — oo electrons are in 
a Landau level, and consequently there is no net cur- 
rent in the y-direction. Decreasing k, results in a net 
current in the y-direction, which is positive for the up- 
per branches (2m + 1), but is initially negative for the 
branches (2m). For more negative values of k it increases 
almost linearly with increasing \k\, which becomes the 
first order analytical result v m (k — > oo) = —k, obtained 
by differentiating (Bah and (|9l 



C. Classical picture 

The center of the classical orbit corresponds to a zero 
in the effective potential. The energy spectrum can be 
divided up into three regions which can be classicaly un- 
derstood by the electron orbits drawn in Fig. |l|(a). In 
region (A) the electrons move in closed orbits either in 
the magnetic field on the left or on the right hand side. 
Since its cyclotron radius is smaller than the distance to 
the magnetic field discontinuity, they feel a homogeneous 
magnetic field. There is no net velocity. In region (B) the 



cyclotron radius intersects the magnetic field discontinu- 
ity slightly, i.e. in such a way that the moving electron 
and the center of its orbit are on the same side. The 
electron is nevertheless able to penetrate in the opposite 
magnetic field region, which results in a (rather small) 
propagation in the y direction v y > 0. For k = the 
center of the orbit is exactly on the edge between the 
two opposite magnetic fields. In region (C) the center is 
located in the opposite magnetic field region of which the 
electron is moving in, resulting in a faster propagation of 
the electron in the y-direction. These different regions 
are also indicated in Fig. ||(a). 

We can also make a quantitative classical study of 
the velocity, starting from the quantum mechanical en- 
ergy spectrum. Since in classical mechanics there is no 
quantization, we make use of the obtained quantum en- 
ergy spectrum in order to find the classical energy and 
thus the radius of the cyclotron orbit. Classically, the 
energy is contained in the circular velocity v v through 
E (k) = v^(k)/2. For any given quantum mechani- 
cal _E(fc)-value we obtain classically the circular velocity 
v<p(k) — y/2E(k). Now if we consider xq — ±k to be the 
center of the electron orbit, we can calculate for every 
k- value the classical velocity v y (k), since we also know 
v v (k) and the cyclotron radius R(k) = v ip (k). Using geo- 
metric considerations, we obtain the following relation 



v y (k) = v v (k)\Jl - [k / v v {k)} 2 / &iccos[k / v^k)], (10) 

which is shown in Fig. ||(a) by the dotted curves. Com- 
paring this with its quantum mechanical counterpart, we 
notice that for k < values good agreement is found, but 
for k > there is a large discrepancy. Moreover, one can 
see that negative velocities cannot exist classically. 

The critical k- value, k*, for which no classical propa- 
gating states can exist, i.e. the electron describes just 
a circular orbit in a homogeneous magnetic field, has 
to be equal to the cyclotron radius k* = R(k*) = 
v v (k*) — ±y^2E (k*), which leads to the boundary 



drawn in Fig. |2j(a) (dotted parabola). 



IV. WITH A BACKGROUND MAGNETIC FIELD 

With a background magnetic field three different con- 
figurations: a) < B < Bo, b) B = Bo, and c) Bo < B 
must be considered. In the following we will study the 
snake orbits in these configurations. 

A. < B < B 

Applying a background magnetic field < B < Bo, 
results in a situation which is very similar to the pre- 
vious one. Again the two magnetic fields have opposite 
sign, but in this case they also have a different strength, 
i.e. B l = —B r /p. Again we can calculate analytically 
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the correction to the energy in the limit k — > oo. For an 
electron on the right hand side in the m th Landau state 
of a magnetic field with strength B r = Bq, the devia- 
tion from the Landau energy due to the presence of the 
other parabola in the effective potential is given by the 
following matrix clement, which to second-order reads 



E m {k 



oo 



\(R\ H \R) (k)\ 



>.!, 



u 2m- 



For an electron on the left hand side, i.e. in the smaller 
magnetic field Br = —Bo/p region, in the m th Landau 
level, this results in 



the left coincides with the first on the right. The classical 
picture for the three different regions corresponds to the 
one drawn in Fig. 0(b), and is also similar to the previous 
case, except for the different cyclotron radii. With this 
picture in mind, one can again calculate the classical ve- 
locity, which turns out to be identical to Eq. ([n]) . From 
Fig. 0(b) we notice that again we obtain good agreement 
for k > 0, but for k < 0, there is a large discrepancy. The 
negative velocity can also in this case not be explained 
H) classically. 

The critical k- value k* = \J2E (k*) for which snake or- 
bits are classicaly possible are indicated by the parabola 
in Fig. 0(b). 



£ m (fc^co) = \(L\H\L) (k)\ 



B. B = B[ 



i) 



im- 2 



'2m- 



Also in this case the energy is smaller then the corre- 
sponding Landau energy. The downward energy shift de- 
creases for increasing p. 

If p is an integer, Landau states on the left and right 
hand side, respectively with quantum number p ■ to and 
to, coincide for k — > oo. As a consequence these states 
have an overlap, which reads to first order 



(L\H\R) =(-!)' 



2 m(p+l)/2 



f 



1/2 



(pm)\m[iT / 

xpP m/2 e -k 2 (l+p)/2 k rn(p+l) + l _ (13) 



One can see that for decreasing magnetic field, i.e. in- 
creasing p, this function decreases because of the expo- 
nential factor. The electron wavefunction in the lower 
magnetic field region is extended over a larger region, and 
further away from the other (center at kp). The overlap 
therefore decreases with increasing p. As a result of this, 
the energy for k — > oo and p > 1 is given by (R\ H \ R) 
and (L\H\L). 

For p = 1, we obtain the previous result, but for in- 
creasing p, the second order term in Eqs. (|ll|) and Jl^ ) 
becomes more important than Eq. (p"3|), because of the 
exponential factor. The splitting is lifted, and the main 
contribution to the negative velocity for k — > oo arises 
from Eq. ([l3|) due to the finite extend of the wavefunc- 
tion in the other parabola. 

As an example we studied numerically the case when 
a background magnetic field B = Bq/2 is applied, i.e. 
B l = —Bq/2 and B r — 3Bq/2. As one can see in 
Fig. 0(b) , this results in two parabolas with different min- 
ima and confinement strength. The resulting spectrum 
(see Fig. 0(b)) is very similar to the one of the previ- 
ous case, but unlike the previous symmetrical case, not 
all states are twofold degenerate for k — > oo. We now ob- 
tain two different sets of Landau states, corresponding to 
electrons moving in different magnetic field regions with 
different strength. In this case the second Landau level on 



(12) when a background magnetic field B = Bq is applied, 
we obtain the magnetic barrier studied in Ref. where 
the magnetic field is different from zero only in the re- 
gion x > 0, i.e. B l = and B r = 2B . From Fig. 0(c) 
one can see that in this case the potential is made up of 
only one parabola and on the left side it is a constant 
k 2 /2. The energy spectrum and corresponding velocities 
for this particular case are shown in Fig. 0(c) and 0(c), 
respectively. We notice that for k — > oo we again obtain 
Landau states, which correspond to bound states on the 
right hand side of the magnetic edge. Consistent, as be- 
ing a limiting case of the former magnetic field states, 
i.e. p = oo, the energy decreases with decreasing k and 
there is no splitting of the energy levels. Thus now we 
only have states which propagate with negative velocity 
to which we cannot assign a classical interpretation. 

Also in this case we can divide up the spectrum into 
three regions: (A) the electrons move in closed orbits 
in the magnetic field region on the right hand side, (B) 
electrons are free, propagate forward and are reflected on 
the barrier and (C) electrons are free, propagate back- 
ward and are reflected on the magnetic edge. Notice that 
for a free electron, the energy is larger than fc 2 /2, since 
now the electron also propagates in the rr-direction and 
consequently has an additional kinetic energy k 2 /2. 

Classically, propagating states in the magnetic field 
region do not exist, only Landau states do. The bound- 
ary where these classical trajectories are possible is again 
given by k* = y/2E(k*). 



C. B > B 

By applying a background magnetic field with strength 
larger than B > Bq, we arrive at the situation where < 
B l < B r . The magnetic fields on the left and the right 
hand side have the same sign, but a different strength, 
i.e. B l = B r /p. 

To obtain the energy in the limits k — > ±oo, we again 
can approximate the wavefunction as being in a Landau 
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state in the corresponding magnetic field. We found 
E{k -400)= (R\ H \R) (k) 

r 

= m+- 

-^Lfl + ^k 2 ™-^-* 2 , (14) 
toIVtt V P ) 

for an electron on the right hand side in the m th Landau 
state of a magnetic field with strength B r = Bq. For an 
electron on the left hand side, in the smaller magnetic 
field B l = Bq/p in the m th Landau level, we have 

E(k^> -oo) = (L\H\L) (k) 

= - m + - + —r-= 1 + - )k 2m l e k , 
P \_ 2 J mi^/TT \ p J 

which results in a negative velocity. 

The energy spectrum and the velocity of these eigen- 
states for the case when B — 3Bq/2, i.e. B l = Bq/2, 
B r = 5B /2, are plotted respectively in Fig. ||(d) and 
||(d) . The center of the orbit is situated on the right side 
for k > 0, for k < it is on the left side. For k — > ±oo, 
the electrons move in a homogeneous magnetic field (on 
the left (k — > — oo) or right (k — > +oo) hand side of 
x = 0), and thus u a = 0. 

From Fig. |l|(d) one notices that there is only one min- 
imum in the effective potential because the minima of 
both parabolas are now situated on the same side. The 
trajectories corresponding with regions (A), (B), and (C) 
are depicted in Fig. 0(d). The trajectories in region (A') 
are similar to those in (A) but now for a magnetic field 
on the left hand side, i.e. with smaller strength. 

Geometrical considerations yield the following classical 
velocity 

v y (k) = 2v v {k)^l-[k/v v {k)] 2 

x {B l arccos [—k/v v (k)] 

+B r arccos[k/v v (k)]}-\ (16) 

which is plotted in Fig. ||(d) as dotted curves together 
with the quantum mechanical group velocity. One can 
see that, in contrast to the previous cases, the negative 
velocity can be understood as classical snake orbits, but 
these snake orbits all run in the same y-direction and now 
there are no states with v y > 0. 

Notice that the quantum mechanical velocity exhibits 
a small oscillatory behaviour on top of a uniform pro- 
file. These whigglcs can be understood from the electron 
distribution over the two parabolas (see Fig. [|). With 
increasing k, the electron distribution is shifted from the 
left parabola to the right one. Due to the wavelike char- 
acter of this distribution, the probability for an electron 
to be in the right parabola (integrated solid region in 
the inset of Fig. ||) exhibits whiggles as function of k, 
with n maxima as shown in Fig. 0. Energetically it is 



favourable for an electron state to have as much as pos- 
sible electron probability in the lower potential region. 
Consequently, when the electron probability in the lower 
potential region attains a maximum, a maximum down- 
ward energy shift will be introduced on top of the overall 
energy change, and this will result in a maximum in the 
group velocity. 

V. NEGATIVE VELOCITY STATE 

Formally, the existence of the quantum mechanical 
negative velocity state can be attributed to the fact that 
shifting two one dimensional potential wells towards each 
other results in a significant rearrangement of the energy 
levels in the composite potential well. Because the com- 
posed well is broader, some states, p.e. the ground state, 
have an energy which is lower than in each of the indi- 
vidual narrower wells. In this particular case, the wave 
vector k measures the distance of the two wells to each 
other, and consequently this energy decrease results in 
a negative group velocity —dE/dk. In this section we 
focus on these negative group velocity states. 

Since the negative velocity states are present for any 
background magnetic field B, but can only be understood 
classically in the situation B > Bq, we will investigate 
the group velocity for a fixed k value with varying back- 
ground magnetic field. In Fig. || the spectrum is plotted 
as function of the applied background field B. We have 
chosen k — 1.5 because in this case a large negative ve- 
locity is obtained for the lowest level when B = 0. Notice 
that for B < Bq: (1) almost all levels decrease in en- 
ergy with increasing background field; (2) there is an 
anti-crossing for E/E = 0.4 + 0.4585/£ (dotted line). 
This anti-crossing occurs when B/Bq = n/(n+ 1), with 
n the Landau level index. For this condition some of 
the Landau levels are degenerate in the limit k — > oo 
(see Fig. ||(b) for the case n = 0); and (3) for B — » B a 
the separation between the levels decreases to zero and a 
continuous spectrum is obtained with a separate discrete 
level at the anti-crossing line. The continuous spectrum 
for B = Bo results from the scattered states in the po- 
tential of Fig. ||(c), while the discrete state is the bound 
state in this potential. For the considered A-value, i.e. 
k = 1.5, only one bound state is found for B = Bq. 

The corresponding group velocity v y = —dE/dk is 
shown in Fig. ^|. Notice that the maximum negative ve- 
locity is obtained near the anti-crossings in the energy 
spectrum (Fig. ||). Near B / Bq = nl (n+ 1) the splitting 
in the energy spectrum (see Fig. ||) is largest and as a 
consequence one of the levels is pushed strongly down in 
energy and consequently v y becomes strongly negative. 
Notice that: (1) every level has some B/B region at 
which v y < 0, and (2) for B / Bq — > 1 the velocity v y — > 0, 
while (3) the enveloppe of {v y ) min in Fig. ^| reaches for 
B /B the v y < value of the B = B state. For B > B 
we have v v < for all states. 
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Using expressions (|o|) and (|l6|), we can also calculate 
the classical velocity corresponding to the energy spec- 
trum in Fig. M. This is shown in Fig. 0. We notice that 
for B/Bq < 1, the classical velocity has a similar be- 
haviour as the quantum mechanical one, except for the 
anti-crossings and the lack of negative velocities, which 
do not have a classical analogon. These negative veloci- 
ties appear suddenly for B/Bq > 1 and exhibit more or 
less the same behaviour. 

As was already apparent from the above study a neces- 
sary condition for the existence of the non-classical edge 
states is the presence of two local minima in the effec- 
tive potential. In the limiting case B = Bq the second 
minima is the limiting case of a flat region in V4 (x) for 
x < 0. But not all these states have a negative velocity. 
How can we classify them? 

From Fig. ||(a,b) one notices that initially (for rather 
small B values) the parabola E = k 2 /2 separates the 
region where only states with positive group velocity ex- 
ist, from the region where also negative velocity states 
are present. This is due to the fact that the value of 
this parabola equals the barrier height between the two 
parabolic wells for the corresponding fc-value. When the 
energy exceeds this barrier, the shape of the wavefunction 
is not determined anymore by the separate parabolas, but 
by the overall composite well width. For decreasing k the 
well is squeezed, and thus all the energy levels are pushed 
upwards, resulting in a positive group velocity. Although 
this is not an exact rule which cannot be extended rig- 
orously throughout the B < B regime, it nevertheless 
provides insight into the fc-values (or B values) for which 
these negative velocities states arise. 

Inspection of the wavefunctions shows that there is 
a feature which marks the negative velocity states, and 
which relates indirectly to the presence of the different 
potential wells. It turns out that if the wavefunction or 
its first derivative exhibits a dip at some x which satis- 
fies <j)(x)<j)"(x) > and (j)'(x) = 0, or <p' (x)<p"'(x) > 
and 4>"{x) = and the condition that <fi(x) ^ 0, then the 
state has a non-classical negative velocity. This is true 
for every k value, as long as B < Bo. This is illustrated 
in Fig. ^ where we plot the wavefunction for k = 1.5, 
n = 2, with background magnetic field B/Bo — 0.68 and 
0.73. The above dip in the wavefunction or its deriva- 
tive (indicated by the dashed circle in Fig. ||) is a result 
of the different potential wells, which have their separate 
influence on the shape of the wavefunction, and therefore 
hamper the matching. The difference in <p{x) being zero 
(or not), can be interpreted as a generalization of match- 
ing the individual states in an asymmetric (symmetric 
way) when the Landau states are degenerate at k — ► oo. 

VI. TIME DEPENDENT CLASSICAL 
INTERPRETATION 

One can make different attempts to link a classical pic- 
ture to quantum mechanics. Often the comparison starts 



with a schematic classical picture which is then supported 
by comparing the quantum mechanical probability den- 
sity with the classical one, obtained through calcula- 
tion of the classical electron trajectory solving Newtons 
equation. For a ID problem one can also verify the clas- 
sical fflptipiT-by inspection of the velocity parallel to the 
edgeJi^llS'EaEl For a cilindrical symmetric problem, the 
classical electron motion can be inferred from the mag- 
netic moment J3ij_the circular current distribution of the 
electron state.LTLj In this paper, a quantitative compar- 
ison was made by use of a quantum mechanical velocity 
parallel to the ID magnetic field discontinuity. In the 
following, we will try a different approach where we will 
construct time dependent states, and in doing so we will 
introduce another feature, i.e. the oscillation frequency 
perpendicular to the magnetic edge. 

A. B = 

We already mentioned before that the solutions for 
this kind of problem are the parabolic cylinder functions 
4>{x) — D E _i [\/2 (x — fc)], matched in such a way that 
we have symmetric and anti-symmetric wavefunctions as 
is shown for k = 2 in Fig. 0(a). At k = oo, these symmet- 
ric and antisymmetric states are twofold degenerate (see 
the two wavefunctions corresponding to the solid square 
in Fig. |l|(a)). Due to this degeneracy any linear combi- 
nation of these states is also an eigenstate. If we take the 
following linear combination 

|m+)=^(|2m) + |2m+l)), 

|m_) = ^(|2m)-|2m + l», [U) 

we arrive at the well known Landau states, i.e. wavefunc- 
tions of electrons located in two different homogeneous 
magnetic field profiles. One electron is moving clockwise, 
while the other is moving counterclockwise. For decreas- 
ing k this degeneracy is lifted. Although taking linear 
combinations of states with a different energy yields a 
time dependent solution, we will extrapolate this picture 
towards all the other states. We choose a new orthonor- 
mal but time dependent basis: 

|m+) = (e' iB2 '"*|2m) + e lE ^ +lt \2m + 1)) /V2 

|m_) = (e iE2 ™ 4 |2m) - e iEam + lt \2m + 1>) /y/2 (18) 

with E m+ = E m _ = (E 2m + E 2m+ i)/2. The re- 
sulting energy spectrum is shown in Fig. ||(a) by the 
dashed curves. The corresponding velocities are plotted 
m Fig. |(a). For every branch there are two states |to+) 
and |to_). 

With these new quantum states much better agree- 
ment is obtained with the corresponding classical re- 
sults (dotted curves in Fig. ^(a)). Because of the ad- 
dition of the two eigenstates the negative velocity almost 
disappeared. Only the lowering of the energy, as was men- 
tioned in the limiting case (i.e. k — > oo), results in a small 



G 



negative velocity, which can't be understood even in this 
picture. Also the boundary which indicates when classical 
states propagate is in much better agreement now. 

Since we now have time dependent states, we can cal- 
culate a new feature: the oscillation frequency oj x in the 
x-direction. The time dependent probability densities of 
the \m+) and |m_) states have the following form: 

\{m+\x) (t)| 2 = i(|(2m|x)| 2 + |(2m+l| a; )| 2 

+2cos[oj x t] (2m\x) (2m + l|x», 
= | (m_|jc) (t + ir/uj x )\ 2 , (19) 

where uj x ,m — (-E^m+i — E2 m )/h is the quantum mechan- 
ical oscillator frequency in the x-direction. 

Classically we can calculate this frequency io x again, 
using simple geometrical considerations, which results in 



arccos(k/v ip (k)) ^ arccos(— k/v v (k)) 



B> 



B' 



ui x (k) 



2 arccos(— k/v v {k)) ' 



(20) 



Both results are plotted in Fig. g(b), and we obtain 
reasonably good agreement between the quantum (solid 
curves) and classical (dotted curves) results. Notice that 
for \k\ > k* , classically ui x = 0, which means that the 
electron does not oscillate between the two different mag- 
netic field regions (i.e. it is not a snake orbit state), but 
it oscillates in a homogeneous magnetic field and conse- 
quently we obtain the time independent eigenstates cor- 
responding to the Landau levels. 

Of course this approach is only useful if proper lin- 
ear combinations are possible. Unfortunately this is not 
the case when a background magnetic field B < Bq is 
applied. 



B. B > B 



The above approach is also fruitful in the case when 



< B l < B r 



We can again add adjacent levels, 



two by two, similar as described before. We can re- 
peat exactly as was done before, and we arrive again 
at the time dependent states of Eq. (|l8|). The energy 
spectrum of these states when B = 3Bq/2 is shown in 
Fig. ||(d), by dashed curves. From Fig. [n](a), we no- 
tice that the classical velocity is in better agreement 
then before, since the amplitude of the whiggles is low- 
ered, due to the summation. The quantum mechanical 
oscillation frequency in the x-direction is again given by 
(E 2m+1 - E 2m ) and plotted in Fig. ^(b). We 
notice that since there are no degenerate states, we al- 
ways have oscillating electrons. For k — > oo the electron 
oscillates with frequency lu x = 2.5lu c , while for k — > — oo 
the electron oscillates with frequency lu x — 0.5lu Ci i.e. 
the electrons circle around in their seperate homogeneous 
magnetic fields. This can also be seen from the classical 
oscillation frequency in the x-direction, which is given by 



(21) 



Notice that also here whiggles in v y are present (see 
Fig. |lO](b)) which are not present in the classical 
results. It is clear that proper linear combinations can 
always be made, as long as B > Bq. 



VII. CONCLUSIONS 

We studied the electron states near discontinuities in 
the magnetic field. Different ID magnetic field profiles, 
i.e. steps, were considered. The quantum mechanical 
energy spectrum was obtained and the group velocity of 
the states was calculated. Their corresponding classical 
orbits were found and the propagating states which are 
located at the magnetic field discontinuity correspond to 
snake orbits. Quantum mechanical magnetic edge states 
were found which move along the magnetic field step in 
opposite direction to the classical snake orbits and which 
cannot be understood classically. We were able to con- 
struct non stationary quantum mechanical states which 
closely approximate the classical solution for the sym- 
metrical case B l = —B r and for the more general case 
B r > B l > B a . 
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FIG. 1. The effective potential (left part of the figure) for 
k = — 2 (dotted curve) and k = 2 (solid curve) for the different 
magnetic field profiles indicated at the top of Fig. ^. The wave- 
functions corresponding to the states in Fig. Pj, as indicated 
by the corresponding symbols, are also shown. Schematic rep- 
resentation of the classical electron trajectories (right part of 
the figure) corresponding to the different regions indicated in 
Fig. ft 



FIG. 2. The energy spectra for different magnetic field pro- 
files (B l ,B r ) shown on top of the figure: (a) (—1,1), (b) 
(-0.5,1.5), (c) (0,2), (d) (0.5,2.5). The dotted curves mark 
the different classical regions. The dashed curves in (a) and 
(b) are the average energy of two adjacent levels. The sym- 
bols correspond to the electron wavefunctions plotted in Fig. 





FIG. 3. The group velocity in the {/-direction (parallel to 
the edge) of the 6 lowest branches of Fig. ^. The dotted curves 
correspond to the velocity calculated for the classical snake 
orbits. 



FIG. 4. The velocity for the n = 4 state of Fig. |(d), and 
the electron probability for the electron to be in the right part 
of the parabola as function of k. The inset shows the effective 
potential and the electron probability at the local maximum 
kl B = 3.4. 



FIG. 5. The energy spectrum at Mb = 1-5 for the 40 lowest 
states as function of the applied background field B/Bq. 



FIG. 6. The group velocity v y as function of B/Bo, corre- 
sponding to the electron states of Fig. ^. 



FIG. 7. The classical velocity in the y-direction as function 
of B/Bo, corresponding to the electron states of Fig. ^| 



FIG. 8. The wavefunctions for Mb = 1-5, n = 2 for two 
negative velocity states (see Fig. [J) with background magnetic 
field B/Bo = 0.68; 0.73. 



FIG. 9. (a) The quantum mechanical velocity in the 
{/-direction (solid curves) of the time dependent states for 
the case B l — —B r = — 1 and the velocity obtained classi- 
cally (dotted curves) . (b) The same as in (a) but now for the 
oscillation frequency in the x-direction. 



FIG. 10. The same as in Fig. |j, but now for the magnetic 
field profile with (B l ,B r ) = (0.5,2.5). 
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